International Journal of Trend in Scientific Research and Development (IJTSRD) 

Volume 4 Issue 3, April 2020 Available Online: www.iitsrd.com e-ISSN: 2456 - 6470 

» ♦ 

Comparative Assessments of Contrast 
Enhancement Techniques in Brain MRI Images 

Anchal Sharma, Mr. Mukesh Kumar Saini 

Department of Electronics & Communication Engineering, 

Sobhasaria Group of Institutions, Sikar, Rajasthan, India 


CLAHE 


I. INTRODUCTION 

As per the latest statistics provided by the WHO, fatal 
injuries by cancer are 8.8 million in worldwide. In the 
literature, brain tumors are classified as malignant cells that 
develop within the brain. These cancer cells develop into a 
mass of cancerous tissues that impede with brain abilities 
comprising of motor function, sensation, memory, and 
various daily body capabilities [1]. 

Many cancerous cells are referred as malignant tumors and 
people formed of explicitly non-cancerous cells commonly 
adverted as benign tumors. Further, there are mainly two 
main types of brain tumors generally referred as primary 
and secondary. Tumors or cancer cells that often emerge 
from the brain tissue are called primary brain tumors, 
whereas tumors that either spread from different brain 
organs are known as secondary or metastatic brain tumors. 

It is possible to remove benign brain tumor, which consists 
of cancer cells. Commonly, benign brain tumors have clear 
margin or edge, generally not expanding to other parts of the 
body. However, benign tumors however, persuade grave 
health challenges. Benign brain tumors are of two type 
Grades I and II [2, 3, 4]. The other form of tumor known as 
malignant brain tumor consists of cancer cells generally 
referred as brain cancer are likely to grow rapidly, and can 
influence ordinary brain tissues near the area. Such sort of 
tumor can be life threatening. Grade III and IV are the grades 
assigned malignant brain tumors [2] [3] [4]. 
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If detected early, many brain tumors are less hazardous and 
almost all are cheaper to take care of. Apparently, we must 
concentrate our resources and address this issue soon as 
possible. MRI image is a user friendly and extensively 
utilized imaging modality for early detection and qualitative 
diagnosis of brain diseases because of its capability to 
outlook numerous human soft tissues/organs with few 
adverse reactions [1]. For image analysis and interpretation, 
two common enhancement techniques were applied, which 
include spatial filtering and shock filtering are evaluated by 
quantifying the image feature through the calculation of the 
MSE and PSNR of images [3, 4, 5]. Once the enhancement is 
done, it can go for further step of segmentation, which will 
be applied with various approaches like thresholding, and 
region based segmentation. 

II. CONTRAST ENHANCEMENT 

The primary purpose of image enhancement is to refine that 
image in such a way that the resulted effect on image will be 
more feasible to diagnosis than the original image for a 
particular application. This improvement mechanism by 
itself may not improve the intrinsic predictive value of the 
data solely; it merely emphasizes certain features of the 
image [5, 6, 7]. For simulation outcomes using MATLAB for 
different contrast enhancement techniques, it is concealable 
that enhancement is solely application based and is well 
demonstrated. We assess the efficiency of two enhancement 
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techniques in this study, based on two parameters of PSNR 
and MSE. 


The primary aim of enhancing contrast is to bring out 
detailed information obscured in an image. For the analysis 
purpose, various contrast enhancement techniques are used 
including Contrast Stretching, Shock Filter, Histogram 
Equalization, Contrast Limited Adaptive Histogram 
Equalization CLAHE. These techniques are applied on the 
three types of brain MRI images, which are normal, benign 
and malignant brain MRI images [2]. Based on two 
evaluation parameters Peak signal to noise ratio [PSNR] and 
Mean square error (MSE] the comparison was made 
between the techniques [4], [11]. It also has been described 
which technique is best suited for brain MRI analysis and 
gives better performance than others give. 


A. Contrast Stretching 

Contrast stretching is a straightforward method for 
enhancing the image contrast by extending the scope of pixel 
intensity values to expand the parameter set necessary. This 
methodology can only apply a linear scaling function to the 
pixel values for images [ 6 ]. 


Through contrast stretching a low-contrast image may be 
converted into a high-contrast image by restoring or 
remapping the gray-level values to the full range of the 
histogram. [13]. It is referred to as dynamic range extension 
in the scope of digital signal processing. This can be 
demonstrated in equation (!]: 



ax, 

P{x-a)+y^, 

Ax-b)+y^, 


Q<x<a 

a<x<b 

b<x<L 



with a minimum principle which keeps the range of the 
filtered image surrounded by the original image range [ 8 ]. 
The dilation and erosion process is expounded pursuant to a 
diminutive time increment dt using a Partial Differential 
Equation (PDE], this creates a sharp discontinuity called a 
borderline shock around two areas of influence and 
ultimately we get a deblurred output. 

The Kramer and Bruckner definition [7] can be describe 
utilizing this subsequent Partial Differential equation (4]: 

= sign{delta{u))\gradienl{u)\ (4] 

Considering a continuous image/(x,y). Then, evolving f 

under the process may generate a class of filtered images. 
The equation (4] can be written as equation (5]: 

Uf = -5'/g/i(Aw)|Aw| (5] 

Where subscript denotes partial derivatives, and 
Vu = {u^,Uyy is the gradient oft/. Let's assume some pixels 

are in the maximum influence zone (negative Laplacian 
i.e. Au = -\-Uyy, is negative which will then an equation ( 6 ] 

is given for the dilation. 

Uf =|Vt/| (6] 

For Laplacian positives, pixels pertain to a minimum zone of 
influence, with Aw < 0, then (5] can be abridged to an erosion 
given by equation (7]. 

M, =-|Vm| (7) 

Therefore the Laplacian's zero-crossings dole out as an edge 
detector. Essentially, the consequence is an 

enhancement/sharpening of image input. 


Where, x and y are input image and the Stretched output 
respectively and a, p and y represents stretching constants, 
acting as factor of multiplier whereas the lower and the 
higher range are represented by a and b while and 
are calculated by equation (2] and (3]: 

ya=aa 

y^= P{b-a)+y^ (- 3 ^ 

The intent of stretching the contrast in the different 
applications is to introduce the image into a scope which is 
more acquainted or normal to the senses, it is therefore also 
called normalization [5]. 

B. Shock Filter 

For deblurring signals and images shock filter is used by 
creating shocks at points of inflection. Shock filters either 
apply erosion or dilation process produces a "shock" 
between two zones of influence, one is for a maximum and 
the other for a minimum signal. 


C. Histogram Equalization 

Typically a histogram represents uniform distribution of 
pixels in a graphical form. The histogram equalization (HE] 
is a widely used image contrast enhancement technique 
because of its simplicity and efficacy [9, 13]. This method 
improves the global image contrast and accommodates 
image intensities to boost contrast by distributing the 
intensity levels that are most widely used [4], thus the 
intensities of the histogram are better distributed. Thus, the 
areas with lower local contrast reach a greater contrast [3], 
[ 10 ]. 


Considering a discrete gray scale image X={x{i,j)} 
comprising of L discrete gray levels expressed 
as{XQX^X 2 .Over a certain imaged, the 

probability density function P{Xj.) is determined as in 
equation ( 8 ]: 

n [ 8 ] 


The premise is that the step of dilation is utilised near a 
maximum, and the process of erosion is used in the vicinity 
of minimum. The determination about pixel's area of 
persuade (whether maximum or minimum] is constructed 
on the Laplacian basis as the pixel is perceived to be a 
maximum for negative Laplacian in the zone of influence, 
and minimum for positive Laplacian. Shock filters comply 


Where, the numbering of occurrence of gray level X 
represented by^k^ ^ total pixel count in the image. 
Defining the cumulative distribution function (cdf] 
consequent to P{Xj.) is mentioned as the following 
equation (9]: 
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C(x)=£p(xJ 

7=0 

Where x is ^^for k = QX,....L-l andO < )< 1. 


Histogram equalization is a mapping mechanism mapping 
the input image throughout the dynamic range {Xq X^_^ ) by 

using the cdf as a level transformation function. A 
transformation function f{x) based on the cdf is defined as 


in equation (10]: 
f{x)=X,+{X^_i-X,).C{x) 


CIO] 


This is the required Histogram Equalization output image. 


D. Contrast Limited Adaptive Histogram Equalization 
(CLAHE) 

Contrast Limited Adaptive Histogram Equalization (CLAHE] 
is an Adaptive Histogram Equalization generalization which 
is used to avert noise amplification problems [6]. This 
CLAHE algorithm separates the images into contextual 
regions and contributes to the equalization of the histograms 
to each [4]. It evinces the distribution of used gray values 
and makes the features of hidden images more evident. The 
approach comes with three parameters: 

> Block size: This represents the scale of the local region 
that equalizes the histogram around a pixel. This scale 
should be greater than the preservable characteristics. 

> Histogram bins: Histogram bins that are utilized in 
numbers for the process of histogram equalization. The 
number of pixels in a block should be smaller than that. 
This value also limits output quantification when 
processing RGB images of 8 bit gray, or 24 bit. 

> Max slope: It restricts the stretch of contrast in the 
feature to pass strength. High local contrast can result in 
very large values. 

The process takes one added ' clip-level' parameter, which 
varies from 0 to 1. This methodology calculates the 
histogram for each pixel, and after that performs the 
equalization operation of the window or block size. 

111. EVALUATION OF PERFORMANCE MATRICES 

The two metrics are utilized to evaluate the performance of 
different methods of enhancement are considered below [4, 
11 ]. 

1. PSNR - The peak signal-to-noise ratio, abbreviated as 
PSNR, is the ratio of the maximum signal power to the power 
of corrupting noise affecting the fidelity of its representation. 
PSNR can be represented by equation (11]. 


PSNR = l0.log^Q 


f 2 A 
Maxj 

MSB 


( 11 ) 


Here, Maxj is the maximum possible pixel value of the 
image. The specimens are shown using linear PCM with B 
bits per sample, Maxj is 2^ -1. 


2. MSE - The Mean Square error abbreviated as MSE 
provides the cumulative squared error between the original 
image and its noisy approximation. The lower the value of 
MSE, the lower the error. MSE is given below by equation 
( 12 ): 


MSE = — 
mn 


m-\ n-\ 


i=0 >0 


[ 12 ] 


Where x{i,j) is noise free mxn gray scale image and y{i,j) 
is noisy approximation of x{i, j ). 


3. Correlation Coefficient - The correlation coefficient is a 
function of correlation connecting two variables, which 
ranges from -1 to 1. The correlation coefficient will be either 
1 or -1 where the two variables are in perfect linear 
connection. The sign varies depending on how the variables 
are related either positively or negatively. If there is no 
linear relation between variables, the correlation coefficient 
is 0. Here two dissimilar types of coefficients for correlation 
are taken into account; first one is the Pearson product- 
moment correlation coefficient which is more widely used in 
measuring the association between two variables, and the 
other being the Spearman rank correlation coefficient, which 
is based on the rank relationship between variables. Given 
paired measurements (XI, Yl], (X2, Y2], (Xn, Yn], the 
Pearson product moment correlation coefficient 
measurement is specified by equation (13]: 

= , (13) 

IV. RESULTS AND ANALYSIS 
Dataset and Software Implementation 

The data set used for the testing and information on the 
implementation of the program are given in the following 
subsections. Experiments were performed to evaluate 
different commonly used enhancement techniques for 
different type of diseased brain MRI images by the 
comparison we can find the best suited method for 
enhancement of Brain MRI image. 

Dataset 

In the current study, 81 patients dataset constituting of 11 
Benign, 25 Gliomas, 30 Meningioma and 15 Metastases, are 
taken from 512 MR brain tumor slices marked by the 
radiologists using CBAC out of which four images from each 
category is shown in the fig. 1 to fig. 4 respectively. These 
images are collected online available dataset from the 
website radiopedia.org. 



(b) 
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(d) 

Fig.l (a), (b), (c) & (d) are different types of Benign 
Brain MRl images 



(b) 




(d) 

Fig. 2 (a), (b), (c) & (d) are different types of Gliomas 
Brain MRI images 




(d) 

Fig. 3 (a), (b), (c) & (d) are different types of 
Meningioma Brain MRI images 
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Software Implementation 

These proposed methods are implemented in MATLAB 9.0 and are tested on various brain tumor MR images of size 400x400. 
The experiments were performed on PC having InteT'' i3 Processor 3.0 GHz processor with 4 GB RAM. The algorithm takes 3 
min for training the samples. 

The comparison of various contrast enhancement techniques for gray scale images is carried out based on the two parameters 
that are Peak Signal to Noise Ratio (PSNRJ and Mean Square Error (MSEJ. These parameters are being used as objective measures 
for evaluating the performance of the improvement methods applied. As per the evaluation, the result of normal, benign and 
malignant brain MRI images are mentioned in following Table I, Table II and Table III. 


TABLE I OUTCOME OF DISTINCT ENHANCEMENT TECHNIQUES FOR BENIGN BRAIN MRI IMAGE 




Parameters 

Image Category 

Technique Applied 

MSE 

PSNR 

cc 

Benign Brain MRI 
Image 

Contrast Stretching 

2051.840 

15.0358 

0.51651 

Shock Filter 

9253.218 

8.94204 

0.58559 

Histogram Equalization 

234.0376 

26.8647 

0.21827 

CLAHE 

995.266 

18.3353 

0.81775 


TABLE III OUTCOME OF DISTINCT ENHANCEMENT TECHNIQUES FOR GLIOMAS BRAIN MRI IMAGE 




Parameters 

Image Category 

Technique Applied 

MSE 

PSNR 

CC 

Gliomas 

Brain MRI Image 

Contrast Stretching 

1047.2492 

18.1140 

0.41745 

Shock Filter 

4928.6040 

11.4430 

0.34158 

Histogram Equalization 

425.00703 

23.6879 

0.21827 

CLAHE 

1287.6206 

17.2889 

0.82576 


TABLE mil OUTCOME OF DISTINCT ENHANCEMENT TECHNIQUES FOR MENINGIOMA BRAIN MRI IMAGE 




Parameters 

Image Category 

Technique Applied 

MSE 

PSNR 

CC 

Meningioma 
Brain MRI Image 

Contrast Stretching 

862.0564 

19.0362 

0.74360 

Shock Filter 

5672.74 

11.6601 

0.27683 

Histogram Equalization 

1073.816 

18.7656 

0.35151 

CLAHE 

1293.88 

17.6669 

0.85294 


TABLE IVV OUTCOME OF DISTINCT ENHANCEMENT TECHNIQUES FOR METASTASES BRAIN MRI IMAGE 



■ ATI 

Parameters 

Image Category 

Technique Applied 

MSE 

PSNR 

CC 

Metastases 
Brain MRI Image 

Contrast Stretching 

1364.4298 

16.9439 

0.29657 

Shock Filter 

5011.9875 

11.3849 

0.52930 

Histogram Equalization 

881.30001 

20.4975 

0.28053 

CLAHE 

3178.0222 

17.6209 

0.88397 
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V. CONCLUSION 

In this research paper, multiple contrast enhancement 
techniques are applied for brain MRI image analysis and 
comparison was made, which is useful in determining the 
best method for clinical diagnosis. It is clear from the above 
comparison tables that the Histogram Equalization filter 
gives the minimum MSE and the highest PSNR value and best 
Correlation coefficient; therefore, it is the best suited method 
and has delivered better performance than others. 
Consequently, this provides sonologist and radiologist better 
visual perception for the diagnostic purpose of brain MRI 
disease. We will be applying other Biomedical Image 
Processing approaches for better performance in future. 
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